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SUMMARY 

(Turbulence,  heat -transfer,  and  boundary  layer  measurements  were  obtained 
in  a  conical  nozzle  operating  in  air  at  a  nominal  total  temperature  and  pres¬ 
sure  of  960°  R  and  300  psia,  respectively.  Experimental  heat-transfer  coeffi¬ 
cients  were  compared  to  values  determined  by  three  prediction  techniques; 
namely,  (l)  Nusselt  number  correlation,  (2)  a  compressible  boundary  layer 
theory,  and  (3)  an  incompressible  boundary  layer  theory.  Methods  (l)  and  (2) 
yielded  nozzle  heat-transfer  coefficients  that  were  appreciably  higher  than  ex¬ 
perimental  values;  however,  method  (3)  predicted  the  nozzle  heat-transfer  coef¬ 
ficients  very  well  yielding  results  that  were  within  about  2  percent  of  the  ex¬ 
perimental  values  at  the  throat^7 

In  part  of  the  investigation  a  simulated  nuclear  reactor  core  turbulence 
generator  installed  3  inches  upstream  of  the  cylindrical  inlet  failed  to  alter 
the  nozzle  heat-transfer  coefficients.  In  addition,  the  turbulence  generator 
had  essentially  no  effect  on  the  boundary  layer  temperature  profiles  in  the 
nozzle. 


INTRODUCTION 

An  accurate  accounting  of  cooling  requirements  for  nozzles  subjected  to 
the  high-energy  gaseous  environments  usually  associated  with  chemical  and  nu¬ 
clear  rockets  necessitates  the  prediction  of  local  gas -side  heat-transfer  coef¬ 
ficients.  Accuracy  in  the  prediction  of  the  convective  heat-transfer  coeffi¬ 
cients  is  especially  important  for  nuclear  rocket  nozzles  since  the  regenerative 
cooling  capabilities  may  be  marginal. 

A  conventional  method  of  predicting  the  local  heat -transfer  coefficients  in 
rocket  nozzles  consists  of  employing  empirical  Nusselt  number  correlations 
(ref.  l),  based  on  some  reference  condition.  These  Nusselt  number  correlations 
completely  neglect  the  effects  of  pressure  gradient,  variation  in  wall  temper¬ 
ature,  and  boundary  layer  growth.  Experience  has  shown  that  Nusselt -type  cor¬ 
relations  do  not  provide  the  accuracy  required  in  the  prediction  of  nuclear 
rocket  nozzle  heat-transfer  rates.  Therefore,  it  is  desirable  to  use  a  more 


fundamental  approach  to  the  determination  of  gas-side  heat-transfer  coeffi¬ 
cients,  which  involves  a  knowledge  of  the  boundary  layer  characteristics.  For 
example,  the  widely  used  compressible  boundary  layer  theory  discussed  in  refer¬ 
ences  2  and  3  provides  a  means  of  predicting  boundary  layer  development  and 
heat  transfer  for  flows  through  convergent-divergent  nozzles. 

The  boundary  layer  theory  of  reference  3  consists  of  a  simultaneous  solu¬ 
tion  of  the  integral  momentum  and  energy  equations.  Although  a  pressure  gra¬ 
dient  term  is  included  in  the  momentum  equation,  a  modified  flat-plate  theory 
was  used  for  the  skin  friction  coefficient  and  Reynolds  analogy.  In  addition, 
l/7-p ower  velocity  and  temperature  profiles  were  used.  An  incompressible 
boundary  layer  theory  (ref.  4)  based  completely  on  flat-plate  theory  is  also 
applied  to  this  conical  nozzle  flow  study. 

In  using  these  boundary  layer  theories  for  predicting  heat  transfer,  an 
estimate  of  the  initial  boundary  layer  conditions  in  the  nozzle  is  necessary. 

The  problem  of  determining  these  initial  conditions  is  compounded  by  factors 
such  as  turbulence  level  and  boundary  layer  transition  phenomena.  Therefore, 
in  experimental  studies  of  heat  transfer,  these  factors  allude  to  boundary 
layer  and  turbulence  measurements  as  well  as  controlled  inlet  flows. 

In  reference  5,  the  results  of  a  nozzle  heat-transfer  investigation  con¬ 
ducted  in  a  hot-air  facility  show  rather  large  differences  between  measured 
heat-transfer  rates  and  values  predicted  by  (l)  a  Nusselt  correlation  and  (2) 
an  application  of  the  boundary  layer  method  of  reference  2  as  modified  in  ref¬ 
erence  6.  The  discrepancy  was  characterized  by  experimental  heat-transfer  co¬ 
efficients  that  were  lower  in  the  throat  section  and  higher  in  the  convergent 
section  of  the  nozzle  than  predicted  values.  These  air-facility  data  tend  to 
parallel  the  heat-transfer  results  obtained  from  actual  firings  of  a  full-scale 
hydrogen-oxygen  rocket  (ref.  7). 

The  present  investigation  was  conducted  to  explore  further  the  nozzle  heat- 
transfer  problem  with  particular  emphasis  on  the  convergent  part  of  the  nozzle. 
This  part  of  the  nozzle  is  of  interest  because  it  is  a  region  of  flow  acceler¬ 
ation  that  influences  the  boundary  layer  characteristics  and,  subsequently,  the 
heat  transfer  in  the  throat  region  of  the  nozzle.  The  primary  objectives  of 
the  study  were  (l)  to  obtain  heat-transfer  and  boundary  layer  measurements  with¬ 
in  the  cylindrical  inlet  and  convergent  portions  of  the  nozzle  and  to  apply  the 
results  in  the  evaluation  of  various  heat-transfer  prediction  methods  and  (2) 
to  observe  the  effects  of  upstream  turbulence  on  nozzle  heat-transfer  rates  and 
boundary  layer  profiles. 

The  experimental  equipment  comprises  a  cylindrical  inlet  followed  by  a 
conical  nozzle  that  was  instrumented  with  static  pressure  taps  and  plug-type 
heat-flux  meters.  A  boundary  layer  suction  mechanism  was  used  to  remove  the 
upstream  boundary  layer  and  to  produce  a  zero  boundary  layer  thickness  at  the 
cylindrical  approach  section.  Boundary  layer  temperatures  and  pressure  surveys 
were  conducted  approximately  midway  in  the  cylindrical  inlet  and  at  one  station 
in  the  convergent  part  of  the  nozzle.  A  hot  wire  located  in  a  plane  1/4  inch 
upstream  of  the  cylindrical  inlet  was  used  to  determine  the  turbulence  level 
for  an  incoming  air  flow,  which  was  maintained  at  a  nominal  stagnation  pressure 
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of  300  pounds  per  square  inch  absolute  and  a  temperature  of  960°  R. 

SYMBOLS 

A  probe  station  (tables  I  and  IV) 

a  heat-meter  thermocouple  nearest  gas  side  of  nozzle  vail 
B  probe  station  (table  i) 

b  heat -meter  thermocouple 

C  probe  station  (tables  I  and  IV) 

Cp  skin  friction  coefficient 

Cp  specific  heat  of  gas  at  constant  pressure 
V  constant  of  integration 

c  heat-meter  thermocouple  nearest  coolant  side  of  nozzle  vail 

D  diameter 

h  heat-transfer  coefficient 

i  enthalpy 

K  thermal  conductivity  of  Inconel 

k  thermal  conductivity  of  air 

k0  intercept  of  conductivity  of  data  at  0°  F 

L  length 

m  slope  of  conductivity  data 

Nu  Russ el t  number 

P  pressure 

Pr  Prandtl  number 

q  heat  flux 

R  inlet  or  nozzle  radius  at  probe  station 

Rep  Reynolds  number  based  on  nozzle  diameter 


Re^  Reynolds  number  based  on  momentum  thickness 
rc  nozzle-throat  radius  of  curvature 

T  temperature 

t  heat-flux  meter  shaft  temperature 

u  velocity 

u+  dimensionless  velocity 

x  axial  coordinate  measured  from  nozzle  throat 

y  distance  from  wall 

y+  dimensionless  distance  from  wall 

3  angular  position  of  nozzle  instrumentation  (table  i) 
€  area  ratio 

A  temperature  difference  ratio 

p  viscosity 

cp  energy  thickness 

p  density 

G  momentum  thickness 

t  shear  stress 

Subscripts : 

ad  adiabatic 

av  average 

e  free -stream  edge  of  boundary  layer 

i  based  on  enthalpy 

0  total  condition  in  plenum 

ref  reference  enthalpy  condition 

s  static  condition 

T  local  stagnation  condition 
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t  throat 
w  wall  condition 
Superscripts : 

time -averaged  value 

'  pertaining  to  fluctuations  from  time  average 


APPARATUS 

The  nozzle  heat-transfer  facility  is  composed  of  (l)  a  heat  exchanger,  (2) 
a  plenum  and  bypass  flow  mechanism,  (3)  a  water-cooled  test  nozzle,  and  (4)  an 
exhaust  system.  A  schematic  diagram  of  the  facility  is  shown  in  figure  1. 


Heat  Exchanger 

A  jet-engine  combustor  was  used  to  heat  the  coils  of  the  heat  exchanger 
through  which  passed  high-pressure  dry  air.  The  air  was  heated  to  a  stagnation 
temperature  of  960°  R  at  a  pressure  of  300  pounds  per  square  inch  absolute. 


Plenum  and  Bypass  Bleed  Mechanism 

The  air  passed  from  the  heat  exchanger  through  a  diffuser  into  a  14 -inch 


White  gas 

fuel-air 

mixture 


High-pressure  air 


_ V  Atmosphere  exhaust 


^Heat  exchanger 


Bypass  bleed 
flow  control 

Turbulence  generator-^  valve-x 
Plenum,  ' 

14-in.  i.d.-t 


(-Cooling 
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i 
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Figure  1.  -  Schematic  diagram  of  nozzle  heat-transfer  facility. 
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diameter  plenum.  A  30-mesh  screen  was  installed  between  the  diffuser  and 
plenum  to  provide  uniform  flow.  The  plenum  was  a  14 -inch-inside -diameter  by 
6-foot-long  chamber  that  converged  to  a  diameter  of  8.4  inches  at  the  cylin¬ 
drical  inlet  of  the  nozzle.  Pitot  pressure  surveys  in  the  plenum  indicated 
that  the  profile  is  radially  uniform  to  within  0.5  inches  of  the  wall.  The 
velocity  in  the  plenum  was  about  17  feet  per  second  with  the  bleed  flow  ad¬ 
justed  to  give  a  uniform  velocity  profile  across  the  cylindrical  approach  sec¬ 
tion.  This  velocity  was  sufficiently  low  to  yield  static -to-total  temperature 
and  pressure  ratios  of  virtually  1.0. 

A  flow  bypass  manifold  was  used  to  remove  the  boundary  layer  along  the 
plenum  wall.  The  bleed-flow  annulus  height  was  0.9  inch  at  the  cylindrical  in¬ 
let,  a  spacing  which  was  confirmed  experimentally  to  be  adequate  for  removal 
of  the  plenum  boundary  layer.  Bleed  flow  rates  were  adjusted  by  means  of  a 
pneumatically  controlled  valve  which  operated  in  a  critical  flow  mode. 


Cylindrical  Inlet  and  Nozzle  Configuration 

An  uncooled  (adiabatic)  approach  section  with  a  17 -inch  length  by  6. 5 -inch 
inside  diameter  was  used  in  the  investigation.  The  cylindrical  inlet  was  ma¬ 
chined  from  AISI  304  stainless-steel  pipe  to  a  wall  thickness  of  0.25  inch. 

The  exterior  surface  of  the  inlet  was  contoured  to  provide  a  leading-edge  ra¬ 
dius  of  0.03  inch. 


In  part  of  the  investigation,  a  turbulence  generator  (fig.  2)  was  installed 
3  inches  upstream  of  the  cylindrical  inlet  to  measure  the  effects  of  upstream 
turbulence  on  the  nozzle  heat-transfer  coefficients.  The  geometry  of  the  tur¬ 
bulence  generator  was  identical  to  the  simulated  nuclear  reactor  core  turbu- 


Figure  2.  -  Upstream  turbulence  generator  (reactor  core  simulator). 


lence  generator  described  in  refer¬ 
ence  5;  however,  in  that  investiga¬ 
tion  the  generator  was  located  0.85 
inch  upstream  of  the  contraction 
section  of  the  nozzle. 

The  cylindrical  approach  sec¬ 
tion  was  attached  to  a  30°-half- 
angle  convergent,  15°-half -angle  di¬ 
vergent  water-cooled  conical  nozzle. 
Flow  disturbances  resulting  from  the 
sharp  corner  formed  by  the  intersec¬ 
tion  of  the  inlet  and  nozzle  were 
assumed  to  be  negligible.  The  ra¬ 
dius  of  curvature  of  the  nozzle  was 
equal  to  the  throat  diameter  (1.492 
in.).  The  exit  diameter  of  the  noz¬ 
zle  was  7.4  inches  corresponding  to 
a  Mach  number  of  approximately  5. 

The  nozzle  configuration  and  a  list 
of  coordinates  are  shown  in  table  I. 
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The  nozzle  was  machined  from  a  forged  billet  of  AISI  304  stainless  steel. 
A  wall  thickness  of  0.5  inch  was  selected  because  of  its  compatibility  with 
the  size  of  the  heat-flux  meters.  Nozzle  cooling  water  was  directed  axially 
over  the  nozzle  by  means  of  a  lucite  shroud.  The  coolant  flow  velocity  was 
maintained  such  that  it  produced  a  net  temperature  rise  of  only  2°  F. 


Exhaust  System 

The  nozzle  was  connected  to  a  24 -inch-diameter  exhaust  duct  by  a  thin  flex 
ible  wafer  flange,  which  yielded  to  nozzle  growth  and  minimized  radial  heat  con 
duct ion.  The  exhaust  pressure,  nominally  2  pounds  per  square  inch  absolute, 
was  generally  sufficient  to  avoid  nozzle  flow  separation  (noted  occasionally  at 
the  downstream  measuring  station) . 


INSTRUMENTATION 

The  basic  instrumentation  was  concerned  with  the  measurement  of  tempera¬ 
tures,  pressures,  and  turbulence  levels.  Temperatures  along  the  heat-flux  me¬ 
ters  were  used  in  the  determination  of  local  heat-transfer  rates  and  wall  tem¬ 
peratures  for  the  nozzle.  Instrumentation  stations  for  the  nozzle  as  well  as 
the  probing  stations  are  listed  in  table  I.  In  addition  to  the  stations  listed 
in  table  I,  five  static  pressure  and  wall  temperature  stations  were  located 
along  the  cylindrical  inlet.  Wall  temperatures  were  measured  on  the  inner  and 
outer  surfaces  of  the  inlet  at  each  station.  Static  pressures  were  measured  on 
the  inner  surface  except  for  the  upstream  station  at  which  static  pressure  meas 
urements  were  obtained  on  both  sides  of  the  leading  edge.  The  leading-edge 
pressure  difference  was  used  to  establish  the  bypass  flow  rate. 


Pressure 

The  plenum  pressure  and  nozzle  static  pressures  8a,  9,  and  10,  were  meas¬ 
ured  with  Bourdon  tube  gages  having  an  accuracy  of  0.25  percent  of  full  scale. 
The  cylindrical  inlet  and  remainder  of  the  nozzle  static  pressures  were  meas¬ 
ured  by  means  of  manometers  connected  differentially.  Manometer  fluids  consist 
ing  of  mercury,  acetylene  tetrabromide,  and  dibutyl  phthalate  were  selectively 
used  to  obtain  the  maximum  sensitivity  at  each  station.  Manometer  fluid  temper 
ature  corrections  were  applied  to  the  data. 

Boundary  layer  total  pressure  measurements  in  the  cylindrical  inlet  and 
nozzle  were  obtained  with  a  miniaturized  pitot  probe  sketched  in  figure  3. 

Prior  to  shaping  the  pressure  sensing  tube,  it  was  necessary  to  vacuum  anneal 
the  0.028-inch-outside-diameter  by  0. 006 -inch-wall  AISI  304  stainless -steel 
hypodermic  tubing.  The  tubing  was  annealed  to  a  hardness  that  permitted  the 
shaping  and  honing  of  a  sharp  inlet  (lip  thickness  of  approx  zero  in.)  without 
cracking  or  tearing.  The  probe  tip  had  a  rectangular  geometry  with  inside  di¬ 
mensions  of  0.002  inch  by  0.030  inch  wide. 

The  probe  traverse  was  terminated  0.002  inch  from  the  wall  by  means  of  a 
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/-•Razor-edge  thick- 


collar  located  on  the  strut  of  the  probe.  When  electrical  contact  was  made  be¬ 
tween  the  collar  and  the  nozzle  wall,  the  resulting  short  circuit  opened  the 
contacts  of  a  relay  which,  in  turn,  stopped  the  actuator  motor.  The  collar  stop 
was  required  to  avoid  contact  between  the  nozzle  wall  and  annealed  tip.  At¬ 
tempts  to  stop  the  probe  without  a  collar  resulted  in  damage  to  the  razor-edge 
tip. 

The  pitot  probes  were  calibrated  for  velocity -head  recovery  factor  over  a 
simulated  boundary  layer  Mach  number  range.  Variation  of  the  Mach  number  up  to 
a  value  of  0.1  had  virtually  no  effect  on  the  recovery  factor.  Recovery  fac¬ 
tors,  which  varied  from  about  97  to  99  percent,  introduced  negligible  error  in 
the  velocity  profile. 

The  boundary  layer  pressure  probes  were  connected  to  dibutyl  phthalate  U- 
tube  manometers  that  were  referenced  to  the  static  pressure  at  the  same  station. 

In  certain  tests  a  modified  pitot  probe  having  a  tempered  stainless -steel 
tip  was  close  coupled  to  a  pressure  transducer  to  obtain  measurements  in  the 
region  near  the  wall.  The  output  of  the  pressure  transducer  was  monitored  on  an 
X-Y  recorder. 
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Figure  4.  -  Inconel  heat-flux  meter  at  nozzle  throat.  (All  dimensions  in  inches.) 


Temperature  and  Local  Heat  Flux 

Steady-state  measurements  of  gas -side  wall  temperature  and  local  heat- 
transfer  rates  were  obtained  by  an  Inconel  plug-type  heat-flux  meter  of  the  type 
described  in  reference  5.  A  sketch  of  the  heat -flux  meter  is  shown  in  figure  4. 
The  0.125-inch-diameter  shaft  contained  three  Chromel-Alumel  0.003-inch  wire 
thermocouples,  spot  welded  in  an  inert  gas  environment  to  the  shaft.  A  calibra¬ 
tion  of  thermal  conductivity  as  a  function  of  temperature  for  the  Inconel  is 
given  in  reference  5. 

(K  =  1.307  X  10-7  t  (°F)  +  1.704  X  10-4  Btu/(sec) (in) (°F)) 

This  equation  represents  the  calibration  data  within  ±4  percent  for  a  tempera¬ 
ture  range  of  200°  to  1000°  F.  Nineteen  of  these  heat-flux  meters  were  in¬ 
stalled  in  the  nozzle,  each  located  180°  from  the  static  pressure  tap  in  the 
same  plane.  Heat -flux  meter  locations  are  given  in  table  I.  The  cylindrical 
inlet  section  inside  and  outside  wall  temperatures  and  the  total  temperature 
of  the  gas  were  measured  with  Chromel-Alumel  thermocouples.  All  the  aforemen¬ 
tioned  temperatures  were  recorded  on  an  automatic  voltage  digitizer  and  paper 
tape  system. 

Plug-type  heat-flux  meters  can  be  subject  to  large  errors  resulting  from 
conduction  through  the  flange  between  the  nozzle  wall  and  the  plug.  The  plugs 
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used  in  this  investigation  were  installed  with  a  light  push  fit  in  order  to  pro¬ 
vide  a  contacting  interface  between  the  plug  and  the  thin  (0.015  in.)  flange  at 
the  gas -side  wall  of  the  nozzle.  The  effectiveness  of  the  contacting  interface 
in  reducing  the  heat  exchange  between  the  plug  and  the  wall  depends  primarily  on 
the  metallic  contacting  pressure  and  gas  pressure  in  the  interface  region.  Es¬ 
timates  of  the  error  in  measured  heat-transfer  coefficients  resulting  from  the 
conduction  heat  transfer  indicate  that  the  experimental  heat-transfer  coeffi¬ 
cients  should  be  multiplied  by  0.9  with  an  uncertainty  of  ±10  percent. 

Total  temperature  profiles  in  the  boundary  layer  were  obtained  with  a  mini¬ 
aturized  probe  of  the  type  shown  in  figure  5.  A  commercial  Chromel-Alumel 
thermocouple  having  a  wire  diameter  of  0.001  inch  was  housed  in  a  conical  after¬ 
body,  which  was  silver  brazed  to  a  0.125 -inch-outside -diameter  support  strut. 

The  metallic  sheath  surrounding  the  thermocouple  was  electrically  insulated  from 
the  body  of  the  probe  so  that  traversing  motion  could  be  stopped  on  contact  with 
the  nozzle  wall  by  the  same  technique  used  for  the  pressure  probes.  The  thermo¬ 
couple  ball  junction  had  a  diameter  of  0.005  inch;  however,  the  sheath  diameter 
of  0.014  inch  limited  the  thermocouple  to  a  minimum  distance  of  0.007  inch  from 
the  wall. 

The  probe  temperature  recovery  factor  was  assumed  to  be  100  percent  since 
Mach  numbers  were  low  (Mach  number,  <  0.08).  Also,  probe  heat  losses  were  as¬ 
sumed  to  be  negligible. 


The  boundary  layer  temperature  probes  were  referenced  to  the  stagnation 


temperature,  and  the  resulting  differential  temperature  signal  was  monitored  on 
an  X-Y  recorder. 


Turbulence 

Turbulence  intensity  profiles  were  obtained  0.25  inch  upstream  of  the  cy¬ 
lindrical  inlet  by  means  of  a  commercial  hot-film  sensor  and- a  constant  temper¬ 
ature  compensating  anemometer.  In  tests  incorporating  the  turbulence  generator, 
centerline  values  of  turbulence  intensity  were  also  measured  in  the  convergent 
part  of  the  nozzle  (Mach  number,  0.08).  The  sensing  element,  comprised  of  a 
platinum  film  on  glass,  had  a  diameter  of  0.001  inch  and  was  designed  for  tem¬ 
perature  environments  up  to  1200°  F.  The  sensor  was  oriented  perpendicular  to 
the  mean  flow  direction  and,  consequently,  provided  a  measure  of  velocity  fluc¬ 
tuations  parallel  to  the  mean  flow  direction. 

Calibration  of  the  sensor  was  obtained  by  varying  the  mean  velocity  at 
constant  temperature  and  plotting  the  sensor  power  as  a  function  of  the  square 
root  of  mass  flux  density.  Five  points  were  generally  used  to  establish  the 
calibration,  which  was  quite  linear  over  the  range  of  flows  investigated. 

PROCEDURE 

All  tests  were  conducted  at  a  nominal  stagnation  pressure  of  300  pounds 
per  square  inch  absolute  at  a  stagnation  temperature  of  970°  R.  Approximately 

l|-  hours  were  required  for  the  system  to  reach  equilibrium.  When  equilibrium 
was  approached,  the  bypass  bleed  flow  rate  was  established  in  order  to  produce 
a  uniform  inlet  velocity  profile.  The  optimum  bleed  flow  rate  was  determined 
by  a  combination  of  three  methods;  namely,  (l)  observation  of  the  direct  cur¬ 
rent  output  of  the  hot  film  sensor  on  an  X-Y  recorder,  (2)  observation  of  the 
pitot  pressure  on  the  X-Y  recorder,  and  (3)  measurement  of  the  static  pressure 
difference  between  the  inner  and  outer  surfaces  of  the  airfoil-shaped  cylindri¬ 
cal  inlet  lip. 

The  length  of  a  typical  steady-state  run  was  about  2  hours,  as  dictated  by 
the  response  time  of  the  boundary  layer  pressure  probes.  All  data,  with  the 
exception  of  the  boundary  layer  pitot  pressures,  were  obtained  in  a  15  minute 
interval.  Each  temperature  was  recorded  10  times  on  the  automatic  voltage  dig¬ 
itizer  in  order  to  confirm  the  established  steady-state  conditions.  Bourdon 
gages  were  read  and  the  manometers  were  photographed  four  times  in  the  same 
period.  During  the  2 -hour  span  when  boundary  layer  pressures  were  obtained, 
stagnation  pressure  drift  was  limited  to  about  ±0.5  pound  per  square  inch  and 
the  stagnation  temperature  varied  by  ±1.5°  F;  however,  during  the  15-minute 
data  recording  interval  these  tolerances  were  reduced  to  approximately  ±0.25 
pound  per  square  inch  and  ±0.5°  F,  respectively. 

EXPERIMENTAL  RESULTS 
Nozzle  Pressure  Distribution 

Static -to-total  pressure  ratios  for  the  conical  nozzle  are  shown  in  fig¬ 
ure  6  and  table  II.  A  characteristic  deviation  from  one -dimensional  values  can 
be  noted  in  the  neighborhood  of  the  throat. 

Experimental  pressure  ratios  are  in  good  agreement  with  the  data  of  refer - 
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TABLE  II.  -  COMPARISON  OF  EXPERIMENTAL  AND  ONE-DIMENSIONAL 


PRESSURE  RATIOS  FOR  30° -15°  CONICAL  NOZZLE 


Pressure 

tap 

Ratio  of 
axial  distance 
to  throat 
diameter, 
x/Dt 

Area 

ratio, 

e 

Pressure  ratio, 

p/p0 

One- 

dimensional 

Experimental 

Percent 

deviation 

1 

-3.026 

17.436 

0.9992 

0.9999 

0.07 

2 

-2.356 

11.593 

.9983 

.9992 

.09 

3 

-1.684 

6.910 

.9950 

.9963 

.13 

4 

-1.446 

5.541 

.9922 

.9939 

.17 

5 

-1.214 

4.351 

.9874 

.9895 

.21 

6 

-.979 

3.299 

.9779 

.9812 

.34 

7 

-.744 

2.385 

.9568 

.9633 

.67 

8 

-.411 

1.382 

.8546 

.8413 

-1.58 

8a 

-.241 

1.116 

.7385 

.7101 

-4.00 

9 

-.117 

1.027 

.6365 

.5723 

-11.2 

10 

.000 

1.000 

.5283 

.4585 

-15.2 

11 

.087 

1.013 

.4500 

.3731 

-20.6 

12 

.171 

1.054 

.3755 

.2778 

-35.1 

13 

.263 

1,150 

.2886 

.2044 

-41.1 

14 

.425 

1.335 

.2035 

.1822 

-11.7 

15 

.818 

1.873 

.1056 

.1108 

4.65 

16 

1.834 

3.639 

.0346 

.0374 

7.50 

17 

3.665 

8.368 

.0096 

.0092 

-4.35 

18 

5.497 

15.008 

.0040 

.0038 

-5.26 

19 

7.394 

23.992 

.0020 

.0022 

10.0 

-4  -2  0  2  4  6  8 

Ratio  of  axial  distance  to  throat  diameter,  x/Dt 


Figure  6.  -  Pressure  ratio  for  conical  nozzle  in  air.  Total  pressure, 
300  pounds  per  square  inch  absolute;  total  temperature,  967°  R; 
throat  diameter,  1. 492  inches. 


ence  8  and  others  for  conical  nozzles 
of  similar  geometry  (fig.  6).  The 
maximum  deviation  of  -41.1  percent 
from  one -dimensional  values  occurred 
at  an  area  ratio  of  1.15  downstream 
of  the  throat. 

Wall  Temperature  Distribution 

The  Fourier  conduction  equation 
for  one -dimensional  steady-state  heat 
transmission  was  used  in  conjunction 
with  the  measurements  obtained  with 
the  heat  flux  meter  of  figure  4  to 
determine  gas -side  nozzle -wall  tem¬ 
peratures  and  local  heat-transfer 
rates.  The  linear  variation  of  ther¬ 
mal  conductivity  with  temperature, 
determined  by  a  calibration  of  the 
Inconel  plug  material,  was  incorpora¬ 
ted  in  the  equation.  Integration  of 
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TABLE  III.  -  HEAT -FLUX  METER  AM)  NOZZLE  WALL  TEMPERATURES 


[Total  temperature,  958.1  °R;  total  pressure,  299.5  psia. ] 


Heat- 

flux 

meter 

Ratio  of 
axial 
distance 
to  throat 
diameter . 
x/Dt 

Heat -flux  meter 
temperature,  T,  °F 

Wall  temperature,  Tw,  °F 

Wall  temperature,  Tw, 
deviation  from 
average,  percent 

Thermocouple 

Thermocouple 

a  and  b 

a  and  c 

b  and  c 

average 

a  and  b 

a  and  c 

b  and  c 

a 

b 

c 

1 

-3.026 

193.6 

159.0 

125.8 

221.0 

223.1 

222.6 

0.5 

-0.7 

0.2 

2 

-2.356 

210.1 

172.4 

— 

— 

— 

3 

-1.684 

245.5 

197.5 

155.7 

289.3 

277.9 

286.6 

284.6 

1.7 

-2.4 

.7 

4 

-1.446 

260.4 

215.8 

168.1 

302.4 

303.0 

302.6 

302.7 

-.1 

.1 

.0 

5 

-1.214 

273.3 

224.7 

175.9 

319.5 

320.8 

319.8 

320.0 

-.2 

.2 

.0 

6 

-.979 

279.0 

227.7 

172.3 

336.1 

340.4 

337.2 

337.9 

-.5 

.7 

-.2 

7 

-.744 

295.0 

238.2 

177.4 

356.7 

349.9 

354.8 

353.8 

.8 

-1.1 

.3 

8 

-.411 

319.4 

259.0 

198.6 

374.3 

370.0 

373.2 

372.5 

.5 

-.7 

.2 

9 

-.117 

319.9 

257.6 

196.0 

372.1 

368.7 

371.3 

370.7 

.4 

-.6 

.2 

10 

.000 

310.8 

252.6 

192.4 

369.3 

361.3 

367.2 

365.9 

.9 

-1.3 

.4 

11 

.087 

308.7 

250.2 

193.1 

363.6 

354.8 

361.4 

359.9 

1.0 

-1.4 

.4 

12 

.171 

302.7 

265.1 

194.0 

338.1 

383.6 

349.9 

357.2 

-5.4 

7.4 

-2.0 

13 

.263 

293.3 

239.4 

181.4 

340.3 

336.5 

339.3 

338.7 

.5 

-.7 

.2 

14 

.425 

280.6 

232.0 

177.4 

327.2 

329.3 

327.8 

328.1 

-.3 

.4 

-.1 

15 

.818 

258.5 

217.9 

166.9 

302.6 

301.5 

302.3 

302.1 

.2 

-.3 

.1 

16 

1.834 

207.9 

170.2 

132.2  . 

238.7 

236.8 

238.2 

237.9 

.3 

-.4 

.1 

17 

3.665 

122.5 

103.1 

86.1 

140.3 

135.4 

139.1 

138.3 

1.5 

-2.1 

.6 

18 

5.497 

84.7 

73.1 

60.8 

95.3 

96.3 

95.5 

95.7 

-.5 

.6 

-.1 

19 

7.394 

78.4 

67.7 

58.8 

88.4 

84.8 

87.5 

86.9 

1.7 

-2.4 

.7 

the  conduction  equation  yields  a  temperature  distribution  of  the  following  form: 


qy  = 


mt 


k0t  + 


(1) 


The  unknowns  q  and  V  were  determined  by  the  simultaneous  solution  of  two 
equations  by  using  the  shaft  temperature  t  and  corresponding  location  y  at 
any  two  stations  along  the  heat-flux  meter.  Subsequently,  the  wall  temperature 
was  computed  by  setting  y  =  0. 

Measurement  of  three  temperatures  along  the  meter  shaft  provided  three 
possible  combinations  of  variables  for  computing  the  nozzle  wall  temperature. 

If  the  arithmetic  average  of  the  three  wall  temperatures  (in  °F)  agreed  within 
±2  percent  of  the  values  computed  by  each  pair  of  temperatures,  only  the  aver¬ 
age  value  was  used  in  subsequent  calculations  of  heat-transfer  coefficients. 

If  any  of  the  three  wall  temperatures  exceed  ±2  percent  of  the  average  value, 
the  three  temperatures  plus  the  average  value  were  used  in  the  computation  of 
heat-transfer  coefficients.  Shaft  and  wall  temperatures  are  listed  in 
table  III  for  a  typical  run  in  which  the  total  temperature  was  958.1°  R  at  a 
total  pressure  of  299.5  pounds  per  square  inch  absolute.  The  wall  temperature 
in  F,  as  computed  from  each  pair  of  thermocouples,  was  generally  within  about 
±1.0  percent  of  the  arithmetic  average  of  the  three  wall  temperatures. 
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Figure  7.  -  Nozzle  wall  temperature  distribution.  Throat  diameter, 
1. 492  inches. 


TABLE  IV.  -  CENTERLINE  TURBULENCE  LEVELS  WITH 
AND  WITHOUT  THE  UPSTREAM  TURBULENCE  GENERATOR. 


Probe 

station:  A 


[Stagnation  temperature,  TQ,  ~  968°  Rj 


stagnation  pressure,  PQ,  300  psia. ] 


Bleed  flow  optimized 

Turbulence  No  turbulence 

Station 

A 

C 

c 

pu 

' 

0.101 

0.018 

0.028 

pu, 

lb/(ft2)(sec) 

39.3 

121.0 

39.0 

pVK 

lb/(ft2) (sec) 

4.0 

2.2 

1.1 

The  wall  temperature  distribu¬ 
tion  for  the  nozzle  is  shown  in 
figure  7.  Data  are  presented  for 
two  runs  in  which  the  stagnation 
temperatures  differed  by  9.3°  R. 

The  sensitivity  of  the  heat -flux 
meters  to  this  modest  difference  in 
stagnation  temperature  is  reflected 
in  a  slight  shift  in  the  wall  tem¬ 
perature  distribution  curve.  The 
ratio  of  peak  wall  temperature 
(839.4°  R)  to  a  stagnation  tempera¬ 
ture  of  967.4°  R  was  similar  to 
that  of  reference  5  (Tw/T q  ~  0.85) 

for  a  nozzle  of  similar  upstream 
geometry  but  operating  at  a  stag¬ 
nation  temperature  of  1600°  R. 


Turbulence  Measurements 

The  ratio  of  turbulence  inten¬ 
sity  to  the  local  time  average 

value  p^F/pu  at  the  c  enter - 

line  l/4  inch  upstream  of  the  plane 
of  the  cylindrical  inlet  leading 
edge  varied  from  0.101  with  the 
turbulence  generator  to  0.028  with 
no  upstream  turbulence  generator 
(table  IV) .  These  values  were  ob¬ 
tained  for  a  bleed  flow  rate  that 
provided  the  most  uniform  inlet 
velocity  profile.  Further  refer¬ 
ence  to  data  in  which  the  turbu¬ 
lence  generator  was  used  will  be 
confined  to  this  optimized  bleed 
flow  rate. 

Turbulence  intensity  profiles 
in  a  plane  l/4  inch  upstream  of  the 
cylindrical  inlet  are  presented  in 
figure  8.  In  the  absence  of  the 
upstream  turbulence  generator,  the 
turbulence  intensity  profile  was 
essentially  uniform;  for  example, 

p|/ (u* )2/pu  ranged  from  about 
0.028  at  the  centerline  to  0.035 
near  the  wall  of  the  cylindrical 
inlet.  With  the  turbulence  gener- 
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Radius,  r,  in. 


Figure  8.  -  Turbulence  intensity  profiles  at  cylindrical  inlet  for  optimized 
bleed  flow  rates. 


at or  installed,  however,  turbu¬ 
lence  levels  approaching  a  value 
of  0.3  were  obtained  near  the  in¬ 
let  leading  edge. 

It  can  also  be  noted  in  fig¬ 
ure  8  that  the  diameter  of  the 
core  of  uniform  turbulence  level 
was  appreciably  less  than  the  in¬ 
let  diameter.  The  inability  to 
achieve  the  desired  uniform  core 
across  the  cylindrical  inlet  was 
attributed  to  blockage  in  the 
bleed  flow  annulus  resulting  from 
the  turbulence  generator  support. 

A  hot  film  sensor  was  also 
installed  at  station  C  in  the 
nozzle  to  ascertain  the  extent  of 
centerline  turbulence  decay  for 
upstream  turbulence  generation 
(table  IV) .  A  comparison  of 
measurements  at  stations  A  and  C 
indicates  that  the  turbulence  in¬ 
tensity  ratio  p  ^  (uT02 /pTI 
diminishes  greatly  as  the  flow 
passes  through  the  cylindrical 
inlet  and  nozzle  entrance.  Val¬ 
ues  of  p/(u02/ptr  decayed 

from  10.1  to  1.8  percent.  The 
level  of  the  fluctuating  compo¬ 
nents  pj/(u^  decreased  by 
about  a  factor  of  two  between 
stations  A  and  G  (table  IV) . 
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Figure  9.  -  Temperature  profiles  at  convergent  nozzle  probing  station 
with  and  without  turbulence  generator. 


Boundary  Layer  Profiles 

Stagnation  temperature  dis¬ 
tribution.  -  The  boundary  layer 
stagnation  temperature  distribu¬ 
tion  for  the  convergent  section 
of  the  nozzle  (Mach  number, 

0.08)  is  presented  in  figure  9 
in  terms  of  the  difference  ratio 
A  as  a  function  of  the  dimen¬ 
sionless  distance  y/cp.  The 
temperature  difference  ratio  A 
and  the  energy  thickness  <p,  are 
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given  by  the  following  equations: 


T„  T 
T  -  w 

T  -  T 
T,e  w 


-j[  -gk1-  * 


In  the  above  expressions,  it  was  assumed  that  Tip  e  ^  Tq  because  the  Mach  num¬ 
ber  was  low  (0.08)  and  that  the  radius  R  of  the  measuring  station  is  much 
larger  than  9. 

The  temperature  distributions  for  the  convergent  nozzle  probing  station  in 
figure  9  are  presented  for  runs  with  and  without  the  upstream  turbulence  gener¬ 
ator  for  a  nominal  total  temperature  and  pressure  of  968°  R  and  300  pounds  per 
square  inch  absolute,  respectively.  The  energy  thickness  (approx  0.006  in.) 
and  temperature  profiles  were  nearly  identical  for  the  two  cases  indicating 
that  the  upstream  turbulence  generator  had  little  effect  on  the  boundary  layer 
temperature  distribution. 

In  the  uncooled  cylindrical  inlet,  the  departure  from  a  uniform  stagnation 
temperature  distribution  was  essentially  negligible  (A  ~  0)  for  all  runs,  with 
and  without  upstream  turbulence  and,  therefore,  will  not  be  presented. 


Velocity  profiles.  -  The  velocity  head 


measured  in  the  boundary 


Cylindrical  inlet  probing  station 
Cylindrical  inlet  probing  station 
Convergent  nozzle  probing  station 
Convergent  nozzle  probing  station 
Modified  probe  data  at  convergent 
nozzle  probing  station 


20  30  40  ! 


Ratio  of  distance  from  wall  to  momentum  thickness,  y/0 


(a)  Without  turbulence  generator.  Stagnation  temperature, 

967. 4°  R;  stagnation  pressure,  300. 3  pounds  per  square  inch 
absolute. 


(b)  With  turbulence  generator.  Stagnation  temperature,  969. 7°  R; 
stagnation  pressure,  296.45  pounds  per  square  inch  absolute. 


Figure  10.  -  Velocity  profiles  at  cylindrical  inlet  and  convergent  nozzle  probing  stations. 
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layer  was  converted  to  velocity  by  locally  applying  the  incompressible  Bernoulli 
equation.  In  the  cylindrical  inlet  the  density  through  the  boundary  layer  was 
assumed  constant  since  the  wall  was  uncooled;  however,  in  the  convergent  part  of 
the  nozzle  where  heat  transfer  was  present,  it  was  necessary  to  account  for  the 
boundary  layer  density  gradient.  In  figures  10(a)  and  (b),  the  dimensionless 

velocity  head  (P«p  -  Ps)/^  Peuf  an(3-  velocity  u/ue  are  presented  as  functions 

u 

of  the  distance  from  the  -wall  y  nondimens ionali zed  by  the  momentum  thickness 
9 ,  where  6  is  given  by  the  following  equation: 

t)  *  (4) 

the  case  of  no  upstream  turbulence 
the  effects  introduced  by  the  sim- 

In  figure  10(a),  the  nondimens ional  velocity  head  curve  for  the  nozzle  con¬ 
vergence  section  includes  data  obtained  with  a  modified  probe  and  recording  sys¬ 
tem.  The  modified  system,  which  permitted  measurements  closer  to  the  wall,  was 
used  to  generate  the  shape  of  the  velocity  head  curve  in  the  region  near  the 
wall. 

A  velocity  profile  corresponding  to  the  commonly  assumed  l/7-power  law  is 
shown  for  the  convergent  station  in  figure  10(a).  Obviously,  the  experimental 
velocity  profile  was  significantly  different  (more  nearly  approximated  by  a 
l/30 -power  law) . 

The  turbulence  generator  produced  only  a  slight  change  in  the  nozzle  con¬ 
vergence  station  velocity  profile;  however,  a  pronounced  difference  in  profile 
shape  can  be  noted  in  the  cylindrical  inlet  station.  This  latter  difference  in 
shape  is  attributed  to  the  nonuniform  inlet  velocity  profile  produced  by  the 
turbulence  generator  (discussed  previously  in  the  section  Turbulence  Measure¬ 
ments  ) . 

In  determining  the  experimental  momentum  thickness  at  the  convergent  meas¬ 
uring  station,  the  velocity  profile  was  extrapolated  to  zero  in  three  different 
ways.  These  extrapolations  (due  to  lack  of  measurements  at  y  <  0.001  in.)  con¬ 
tributed  an  error  of  approximately  20  percent  in  the  momentum  thickness.  Al¬ 
though  this  error  is  significant,  a  greater  error  in  the  momentum  thickness  re¬ 
sults  from  the  almost  flat  velocity  profile,  which  makes  it  extremely  difficult 
to  determine  the  boundary  layer  termination  point.  The  importance  of  this 
boundary  layer  termination  is  further  increased  by  the  fact  that  the  contribu¬ 
tion  at  this  outer  region  of  the  boundary  layer  to  the  momentum  thickness  is 
very  significant.  The  integrand  of  the  momentum  thickness  through  the  boundary 
layer  gives  a  graphical  visualization  of  this  contribution  (outer  boundary 
layer  region)  to  the  momentum  thickness  (fig.  11 ) .  In  order  to  show  the  greatly 
changed  character  of  the  boundary  layer  in  the  convergent  region  compared  with 
the  cylindrical  inlet  section,  the  integrand  of  the  momentum  thickness  at  the 


The  velocity  profiles  in  figure  10(a)  are  for 
generation,  and  those  of  figure  10(b)  include 
ulated  reactor  core  turbulence  generator. 
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Figure  11.  -  Variation  of  integrand  of  momentum  thickness  across  boundary  layer 
for  cylindrical  inlet  and  convergent  nozzle  probing  stations. 

probing  station  in  the  cylindrical  inlet  section  is  also  plotted  in  figure  11. 
With  these  difficulties  considered  (lack  of  data  near  the  wall  and  inaccurate 
boundary  layer  termination),  approximately  80  to  120  percent  error  in  the  mo¬ 
mentum  thickness  can  be  expected  for  the  convergent  nozzle  probing  station. 

Experimental  values  of  momentum  thickness  are  listed  in  figure  10.  Mo¬ 
mentum  thicknesses  for  the  cylindrical  inlet  were  nearly  an  order  of  magnitude 
greater  than  the  values  for  the  convergence  probing  station.  Cold  flow  tests 
produced  no  significant  change  in  the  momentum,  displacement,  or  boundary  layer 
thickness  at  the  convergent  section  probing  station  compared  with  the  hot  flow 
tests.  Therefore,  the  greatly  reduced  boundary  layer  thickness  and  drastically 
changed  velocity  profile  in  the  convergent  section  compared  with  the  cylindri¬ 
cal  inlet  section  can  be  attributed  primarily  to  acceleration  effects.  Addi¬ 
tion  of  the  turbulence  generator  produced  little  change  in  momentum  thickness 
at  both  stations. 

In  figure  12,  the  shear  velocity  profile  for  the  cylindrical  inlet  station 
is  shown  in  terms  of  u+  and  y+,  which  are  given  in  reference  9  as  follows: 
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Figure  12.  -  Shear  velocity  distribution  at  cylindrical  inlet  probing  sta¬ 
tion.  Stagnation  temperature,  967. 4°  R;  stagnation  pressure,  300.3 
pounds  per  square  inch. 
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The  skin  friction  coefficient  Cj 

was  determined  from  the  experi¬ 
mental  data  in  conjunction  with 
the  von  Karmen  equation  for  the 
"law  of  the  wall"  (ref.  10), 
which  can  he  written  in  the  fol¬ 
lowing  form: 


(7) 


The  experimental  data  in  figure  12  are  in  good  agreement  with  the  von 
Karman  form  of  the  law  of  the  wall  up  to  a  value  of  jr  ~  1000  thus  confirming 
the  presence  of  a  turbulent  boundary  layer  in  the  cylindrical  approach  section. 
In  the  outer  part  of  the  boundary  layer  (y+  >  1000)  a  departure  from  the  law  of 
the  wall  can  be  noted.  This  departure,  referred  to  as  the  "law  of  the  wake" 
(ref.  ll),  is  characteristic  of  turbulent  boundary  layers. 


The  data  of  figure  12  tacitly  indicate  that  boundary  layer  measurements 
could  not  be  obtained  in  the  thin  sublayer  that  terminates  at  y+  ~  26. 


Heat-transfer  coefficient.  -  Experimental  values  of  the  heat -transfer  co¬ 
efficient  based  on  enthalpy  were  determined  from  the  following  equation: 


h,  =  T 


1ad  “  1w 

where  the  adiabatic  enthalpy  ia^  was  given  by 

.1/3  /, 


'ad 


i  +  Pr 
s 


(i0  -  1S) 


(8) 


(9) 


The  Prandtl  number  was  assumed  constant  at  0.7.  The  static  and  total  enthalpies 
is  and  1q,  respectively,  were  computed  by  an  iteration  procedure  by  using  the 

equations  of  reference  12  for  entropy,  .enthalpy,  and  specific  heat. 
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Figure  13.  -  Variation  of  experimental  heaMransfer  coefficient  along  nozzle  and 
comparison  with  theoretical  predictions.  Throat  diameter,  1. 492  inches. 


Experimental  values  of  the  heat-transfer  coefficient  are  presented  in  fig¬ 
ure  13  and  table  V  for  runs  with  and  without  the  reactor  core  turbulence  gener¬ 
ator.  Also  shown  are  heat-transfer  distributions  resulting  from  various  theo¬ 
retical  prediction  techniques;  however,  discussion  of  these  curves  will  be  de¬ 
ferred  until  the  following  section.  The  experimental  heat -transfer  distribu¬ 
tions  of  figure  13  were  duplicated  to  within  about  ±2  percent  in  several  runs; 
however,  for  reasons  of  clarity,  only  the  data  for  two  runs  without  upstream 
turbulence  generation  and  one  run  with  the  simulated  reactor  core  turbulence 
generator  are  presented. 


The  turbulence  generator,  installed  3  inches  upstream  of  the  cylindrical 

inlet,  produced  an  inlet  turbulence  level  p  y (u')2/p  u  of  10.1  percent,  as  com¬ 
pared  with  a  2.8-percent  value  without  the  generator.  Although  a  decay  in  these 
normalized  values  of  turbulence  was  encountered  as  the  flow  passed  through  the 
inlet  and  nozzle  entrance  regions,  the  absolute  turbulence  intensity 

at  station  C  was  about  a  factor  of  two  higher  than  the  upstream  value  with  no 
generator.  As  noted  in  figure  13,  heat-transfer  coefficients  were  essentially 
unaltered  even  with  the  higher  levels  of  turbulence  intensity  (table  IV) . 
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The  maximum  heat-transfer  coefficient,  occurring  slightly  upstream  of  the 
geometric  throat,  was  about  0.0044  pound  per  second  per  square  inch  with  and 
without  the  upstream  turbulence  generator.  It  is  interesting  to  note  that  the 
maximum  heat-transfer  coefficient  in  this  investigation  was  approximately 
80  percent  of  the  value  obtained  in  reference  5  for  a  nozzle  with  the  same  con¬ 
vergence  and  throat  geometry.  Both  nozzles  operated  at  the  same  nominal  total 
pressure  (300  lb/sq.  in.  abs),  but  the  nozzle  of  reference  4  was  tested  at  a 
total  temperature  of  1600°  R  rather  than  960°  R.  The  higher  temperature  re¬ 
sulted  in  lower  values  of  mass  flux  pu.  Therefore,  from  a  strict  mass  flux 
consideration  in  which  h<x(pu)0,8,  the  values  of  heat-transfer  coefficient  in 
this  investigation  would  be  expected  to  slightly  exceed  those  of  reference  5. 

It  is  expected,  however,  that  the  differences  in  the  heat-transfer  coefficient 
may  arise  because  of  the  influence  of  upstream  geometry  and  acceleration  on  the 
boundary  layer  growth  at  the  nozzle  entrance. 


COMPARISON  WITH  THEORETICAL  METHODS 

The  experimental  data  are  compared  with  results  from  the  three  following 
calculation  methods: 

(1)  Nusselt  number  correlation 

(2)  Compressible  boundary  layer  theory  (ref.  3) 

(3)  Incompressible  boundary  layer  theory  (ref.  4) 

Values  of  the  heat-transfer  coefficients  obtained  by  the  above  methods  are 
given  in  table  V. 

Hie  heat-transfer  coefficient  calculated  from  the  Nusselt  number  correla¬ 
tion  and  evaluated  at  the  reference  enthalpy  is  given  by 


Wj^ref 

12  Cp^refD 


(10) 


The  subscript  ref  indicates  the  the  flow  properties  are  evaluated  at  the  ref¬ 
erence  enthalpy  condition  and  the  Nusselt  number  ^uref: 


Nuref  = 


1/3 


(11) 


where  the  Prandtl  number  Rr  is  assumed  constant  at  0.71  and  the  Reynolds  num¬ 
ber  Re^  based  on  the  diameter  of  the  nozzle  at  each  station  and  evaluated 

at  the  reference  enthalpy  is  given  by 


Re 


D,ref 


prefueD 

^ref 


(12) 
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TABLE  V.  -  EXPERIMENTAL  AND  PREDICTED  NOZZLE  HEAT-TRANSFER  COEFFICIENTS 
[Total  temperature,  ~  960.0°  R;  total  pressure,  300.0  psia. ] 


Ratio  of 
axial  distance 
to  throat 
diameter, 
x/Dt 

Heat-transfer  coefficient  based  on  enthalpy, 
h-p  lb/(sec)  (in.2) 

Experiment 

with 

turbulence 

generator 

Experiment 

without 

turbulence 

generator 

Nusselt 

number 

correlation 

Boundary  layer  theory 
(ref.  3) 

Arbitrary 

wall 

temperature 
method 
(ref.  4) 

Cylindrical 

inlet 

starting  point 

Nozzle 

starting 

point 

-3.026 
-2 . 356 
-1.684 
-1.446 
-1.214 

10.40X10"4 

13.19 

17.69 

20.37 

24.02 

9. 63X10 ^ 

10.48X10-4 

10.94 

15.49 

18.23 

22.26 

4. 83X10 “4 
9.42 

15.63 

18.25 

23.92 

16. 81X10  "4 
19.46 

22.35 

21.17X10-4 

27.21 

15. 90x10 ”4 
18.01 

21.07 

-.979 

-.744 

-.411 

-.117 

.000 

28.2IX10"4 

34.93 

42.42 

27.32X10"4 

34.87 

42.95 

43.84 

43.61 

27 .85X10 “4 
36.49 

65.28 

75.75 

72.78 

33. 84X10 “4 
55.84 

68.76 

64.59 

54.17 

36. 79X10 ™4 
51.67 

77.92 

85.01 

75.84 

25. 50X10 "4 
39.11 

38.22 

43.36 

42 . 95 

.087 

.171 

.263 

.425 

.818 

39. 74X10 "4 
36.48 

32.53 

25.65 

42. 17X10 “4 
41.03 

36.63 

32.59 

25.43 

67.74X10’4 

59.01 

49.65 

46.25 

33.48 

47.09X10"*4 

42.92 

38.96 

35.25 

29.79 

68 . 34X10 "4 
56.25 

55.42 

52.09 

40.09 

41. 59X10 “4 
34.83 

34.40 

33.74 

26.75 

1.834 

3.665 

5.497 

7.394 

15 . 11X10  “4 
4.34 

2.88 

14 . 74X10“4 
4.80 

2.77 

2.23 

15.58X10"4 

5.44 

2.68 

1.75 

14. 83X10 

5.75 

2.69 

21. 04X10  ”4 
9.33 

6.46 

14. 65X10 "4 
5.44 

3.15 

In  figure  13  the  Nusselt  number  correlation  based  on  reference  enthalpy 
gives  approximately  a  70 -percent  greater  heat-transfer  coefficient  at  the  noz¬ 
zle  throat  than  that  given  by  the  experimental  data.  Similarly ,  the  compress¬ 
ible  boundary  layer  calculation  (ref.  3)  predicts  a  60-percent  higher  heat- 
transfer  coefficient  near  the  throat. than  that  of  the  experimental  data.  This 
calculation  starts  at  the  probing  station  in  the  cylindrical  inlet  by  using  the 
boundary  layer  thicknesses  measured  there.  Also,  a  similar  calculation,  which 
starts  at  the  convergent  nozzle  probing  station  and  uses  the  measured  boundary 
layer  thicknesses,  predicts  a  93-percent  higher  heat-transfer  coefficient  near 
the  throat  than  that  of  the  experimental  data. 

This  theoretical  boundary  layer  calculation  (ref.  3)  is  based  on  flat 
plate  (zero  acceleration)  heat-transfer  theory.  Therefore,  not  only  the  1/7- 
power  velocity  profile  but  also  the  skin  friction  coefficient  and  Reynolds  anal¬ 
ogy  are  based  on  the  flat  plate  theory.  Comparison  of  this  compressible  bound¬ 
ary  layer  theory  with  the  experimental  data  at  the  convergent  probing  station 
gives  the  following  results: 

(a)  The  experimental  velocity  profile  was  not  represented  by  the  usually 
assumed  l/7-power  law  (fig.  10(a)). 
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(b)  The  experimental  momentum  thickness  is  a  factor  of  4  less  than  the 
calculated  value  for  a  free  stream  Mach  number  of  0.08  based  on  reference  3. 

(c)  The  experimental  displacement  thickness  is  about  a  factor  of  10  less 
than  the  calculated  value  based  on  reference  3. 

This  comparison  indicates  that  the  l/7-power  profile  does  not  describe  a 
nozzle  boundary  layer  distribution;  however,  the  sensitivity  of  the  integral 
calculation  of  heat  transfer  to  the  profile  shape  is  questionable.  In  refer¬ 
ence  13,  boundary  layer  shape  parameters  were  relatively  insensitive  to  pro¬ 
file  powers  of  l/5  to  l/ll;  however  profiles  approaching  a  l/30  power  law,  as 
experienced  in  the  present  study,  may  have  an  appreciable  influence  on  the  heat 
transfer.  Also,  significant  errors  in  the  flat-plate  Reynolds  analogy  as 
applied  to  nozzle  flows  will  greatly  change  the  heat  transfer.  This  is  pri¬ 
marily  caused  by  the  direct  sensitivity  of  the  Reynolds  analogy  and  skin  fric¬ 
tion  coefficient  in  the  integral  calculation  of  the  heat  transfer. 

The  incompressible  boundary  layer  theory  of  reference  4,  which  is  based  on 
incompressible  flat  plate  theory  but  includes  variations  in  wall  temperature, 
predicted  the  heat-transfer  coefficient  very  well.  Heat-transfer  coefficients 
at  the  throat  were  within  about  2  percent  of  the  experimental  results;  however, 
this  excellent  correlation  with  experimental  data  may  be  a  result  of  compen¬ 
sating  errors  in  the  flat  plate  theory  as  applied  to  the  present  nozzle  geom¬ 
etry  and  flow  conditions.  Compensating  errors  are  expected  especially  since 
the  predicted  ratio  of  the  thermal  to  velocity  boundary  layer  thickness  at  the 
convergent  nozzle  probing  station  is  approximately  150  percent  greater  than  the 
measured  value.  Therefore  caution  is  suggested  when  applying  this  heat- 
transfer  calculation  method  to  nozzle  flow. 


SUMMARY  OF  RESULTS 

Experimental  velocity  and  temperature  profiles  in  the  boundary  layer, 
heat -transfer,  and  turbulence  measurements  have  been  presented  for  a  conical 
nozzle  configuration  with  a  cylindrical  approach  section  operating  in  air  at  a 
nominal  total  temperature  and  pressure  of  960°  R  and  300  pounds  per  square  inch 
absolute,  respectively.  An  analysis  of  the  results  of  this  investigation  has 
indicated  the  following: 

1.  Experimental  heat-transfer  coefficients  could  not  be  correlated  by 
either  of  the  following  methods:  (l)  a  Nusselt-type  pipe  flow  correlation  and 
(2)  a  compressible  integral  boundary  layer  prediction  technique.  Predicted 
values  of  the  heat-transfer  coefficient  were  significantly  higher  than  the  ex¬ 
perimental  values  in  the  convergent  and  throat  region  of  the  nozzle. 

Application  of  an  incompressible  boundary  layer  technique  resulted  in  good 
agreement  between  predicted  and  experimental  nozzle  heat-transfer  coefficients. 
Predicted  values  in  the  throat  region  were  within  about  2  percent  of  the  ex¬ 
perimental  results;  however,  the  agreement  may  be  fortuitous  since  the  predic¬ 
ted  ratio  of  the  thermal  to  boundary  layer  thickness  at  the  convergent  nozzle 
probing  station  was  approximately  150  percent  greater  than  the  measured  value. 
Therefore  caution  should  be  observed  with  regard  to  a  generalization  of  these 
results . 
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2.  The  introduction  of  a  simulated  reactor  core  turbulence  generator  up¬ 
stream  of  the  cylindrical  inlet  resulted  in  very  high  levels  of  inlet  turbulence 
(ratio  of  root-mean-square  velocity  to  time-averaged  value  >  0.10)  as  well  as  a 
nonuniform  turbulence  profile;  however,  these  normalized  turbulence  levels  were 
much  lower  in  the  nozzle  convergence  section  (~0.02).  Consideration  of  only 
the  fluctuating  component  of  turbulence  intensity  in  the  convergent  part  of  the 
nozzle  reveals  that  the  value  was  about  two  times  the  upstream  value  with  no 
turbulence  generator.  The  effect  on  nozzle  heat-transfer  coefficients  resulting 
from  this  disturbed  flow  was  negligible. 

3.  The  experimental  momentum  thickness  was  a  factor  of  4  less  than  the 
value  calculated  by  a  compressible  boundary  layer  prediction  technique.  The 
corresponding  displacement  thickness  was  a  factor  of  10  less  than  the  predicted 
value . 

4.  The  greatly  reduced  boundary  layer  thicknesses  and  much  steeper  velocity 
profiles  in  the  convergent  section  compared  with  those  of  the  cylindrical  inlet 
section  can  be  attributed  primarily  to  acceleration  effects. 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  October  11,  1965. 
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